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We study the vortex pinning effect in a Bose-Einstein Condensate in the presence of a rotating 
lattice potential by numerically solving the time-dependent Gross-Pitaevskii equation. We consider a 
triangular lattice potential created by blue-detuned laser beams. By rotating the lattice potential, we 
observe a transition from the Abrikosov vortex lattice to the pinned vortex lattice. We investigate the 
transition of the vortex lattice structure by changing conditions such as angular velocity, strength, 
and lattice constant of a rotating lattice potential. Our simulation results clearly show that the 
lattice potential has a strong vortex pinning effect when the vortex density coincides with the 
density of the pinning points. 

PACS numbers: 



I. INTRODUCTION 

Quantized vortices are one of the most character- 
istic manifestations of superfluidity associated with a 
Bose-Einstein condensate (BEC) in atomic gases 
Formations of triangular Abrikosov vortex lattice in 
Bose condensates have been clearly observed by rotating 
anisotropic trap potentials [2|, 0, [j] . Microscopic mech- 
anism of the vortex lattice formation has been exten- 
sively studied both analytically and numerically using 
the Gross-Pitaevskii equation for the condensate wave- 
function I, i, 0, B i, [lO, [O, [il . 

Another interesting development in ultracold atomic 
gases is creating periodic potential using optical lat- 
tices [l3|. The studies of a BEC in optical lattices have 
found a lot of interesting phenomenon, such as the tran- 
sition from the superfluid to Mott insulator phase [14]. 

More recently, combining the above two systems, i.e., 
rotating BEC and optical lattices has attracted growing 
attention. In particular, the vortex phase diagrams of 
a BEC in a rotating optical lattice potential have been 
theoretically studied, since one expects structural phase 
transition of vortex lattice structures due to vortex pin- 
nms: IT^ , [itI , 0] . In the condensed matter systems 
such as superconductors, the vortex pinning due to im- 
purities and defects in solids have been extensively inves- 
tigated [l^, [^. In the atomic condensates, rotating op- 
tical lattices have been experimentally realized at JILA 
group, making ^use of a laser beam passing through a 
rotating mask [21]. This experiment observed the struc- 
tural phase transition of vortex lattice structures in ro- 
tating Bose condensated gases due to vortex pinning by 
the laser beam. In this system, one can control the pin- 
ning parameters by changing conditions such as angular 
velocity, strength, and lattice constant of rotating opti- 
cal lattices. This allows one to investigate vortex pinning 
effect in a quantitative manner.. 

In this paper, stimulated by the JILA experiment, we 
study the vortex pinning effect by numerically solving the 



time-dependent Gross-Pitaevskii equation. We observe a 
structural phase transition of vortex lattice structures of 
a BEC in a rotating triangular lattice potential created 
by blue-detuned laser beams. 



II. FORMULATION 

We consider a Bose-condensated gas trapped in a har- 
monic potential and a rotating periodic potential with an 
angular velocity Vt. The dynamics of the condensate is 
described by the time-dependent Gross-Pitaevskii equa- 
tion. Assuming a pancake-shaped trap, we use a two- 
dimensional Gross-Pitaevskii equation. In the rotating 
frame, the Gross-Pitaevskii equation is given by (for a 
detailed derivation, see for example Ref. |7|]) 



dt 



2m 



V^ + Kxt(r)+^|V^(r,t)r 



V^(r,t). (1) 



Here, Fext(r) = mujl{x^ + y^)/2 + Mattice(r) describes 
the total external potential, Lz = —ih{xd/dy — yd/dx) 
denotes the z component of the angular momentum op- 
erator, g = ATrfi^as/m is the strength of interaction with 
ttg being the s-wave scattering length, and 7 is the phe - 
nomenological dissipative parameter [22I, [23I, [2^, [25l,l26t . 
The lattice potential created by blue-detuned laser beams 
arranged in the lattice geometry is expressed as 



Via 



exp 



(<t/2)2 



(2) 



where rni,n2 = niai + n2a2 describes lattice pin- 
ning points and Vq describes the strength of the laser 
beam. We consider the triangular lattice geometry 
with two lattice unit vectors given by ai = a(l,0), 
a2 = a(-l/2,73/2), where a is the lattice constant. 
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FIG. 1: Density profile (left) and Structure factor profile 
(right) of a BEG in a triangular lattice potential without ro- 
tation by setting Q = 0. The lattice potential geometry is 
triangular lattice at a/aho — 2.2 and Vq/Z^cjo = 5.0. 



Throughout this pa per, we scale the length and energy 
by ttho = yh/rnuoo and hwo- We set the dimensionless 
interaction strength as C = AirasN/R = 1000, where N 
is the total particle number and R is size of the conden- 
sate along the z-axis direction |7|], and the width of the 
laser beam as cr/aho = 0.65. In our parameter set, the 
healing length is ^/aho = 0.12. We numerically solve the 
Gross-Pitaevskii equation in Eq. ([1]) using Fast-Fourier- 
Transform (FFT) technique. 

We first determine the ground-state condensate wave- 
function without rotation by setting 1] = in Eq. ([T]). 
Starting with this wavefunction as a initial state, we nu- 
merically evolve the Gross-Pitaevskii equation with a ro- 
tating lattice potential with a fixed angular velocity 
until the system relaxes into equilibrium. We investigate 
the vortex lattice structure from the equilibrium conden- 
sate density profile n(r) = |?/^(r)p. In order to quan- 
tify the vortex lattice structure, we calculate the density 
structure factor defined by 



^(k) = J drn{r)e 



-ikr 



(3) 



The structure factor given in Eq. ([3]) contains informa- 
tion about the periodicity of the condensate density. For 
triangular lattice geometry, 6'(K) exhibits periodic peaks 
of the regular hexagonal geometry. In order to distin- 
guish between the Abrikosov lattice and the pinned vor- 
tex lattice, we calculate the peak intensity of the struc- 
ture factor |5'(K)| at the lattice pinning point, namely, 
Ki = 27r/a(l,l/V3), K2 = 27r/a(0, 2/^3). In addition 
to the structure factor, we also calculate the lattice po- 
tential energy defined by 



^lattice = J c^r?/;*(r)l/iattice(r)V^(r), 



(4) 



which can also be used to quantify the transition of the 
vortex structure. As we see below, the vortex pinning is 
signified by a marked decrease of £^iattice- 



III. SIMULATION RESULTS 

In this section, we show our numerical simulation re- 
sults for vortex lattice structures of a BEG in a rotating^ 




FIG. 2: Density profiles (left) and Structure factor pro- 
files (right) in equilibrium state after rotating condensates 
with a fixed angular velocity Q/cJo = 0.70 at a/aho = 2.2. The 
figures show vortex lattice structures for different values of the 
strength of a triangular lattice potential; (a) Vo/huJo = 0.1, 
(b) Vo/hcvo = 0.3, (c) Vo/huo = 0.5. 



triangular lattice potential. Fig. 1 shows the equilibrium 
condensate density profile and the density structure fac- 
tor profile in the presence of a triangular lattice potential 
without rotation. 

Fixing the lattice constant with a/aho = 2.2, we in- 
vestigate the transition of the vortex lattice structure by 
changing the strength Vq for various angular velocities 
ft. Fig. 2 shows the density profiles and the structure 
factor profiles in equilibrium state after rotating conden- 
sates with a fixed angular velocity Q/coq = 0.70. One 
can see that for weak lattice potential Vq/Jiloo = 0.1 
(Fig. 2 (a)), vortices form the Abrikosov lattice, which is 
incommensurate with the triangular lattice potential. By 
slightly increasing the lattice strength as Vo/hcoo = 0.3 
(Fig. 2 (b)), vortices start to being partially pinned by 
the triangular lattice potential. For strong lattice poten- 
tial Vo/huJo = 0.5 (Fig. 2 (c)), all vortices are pinned by 
the triangular lattice potential. We thus observed a tran- 
sition of the vortex lattice structure from the Abrikosov 
lattice to the pinned lattice. 

In Fig. 3, we plot |5'(K)| and £^iattice against the lat- 
tice strength Vq. As shown in Fig. 3 (a), for Q/coq = 
0.70, the lattice potential energy £^iattice decreases grad- 
ually, which indicates that vortices are partially pinned. 
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FIG. 3: Lattice potential energy (a) and peak intensity of the 
structure factor at the lattice pinning point (♦) with a/aho = 
2.2 corresponding to the angular velocity; (a) Q/coo = 0.70, 
(b) Q/ooo = 0.55. Here ^latticeo and |5'o(K)| is the lattice 
energy and the peak intensity of the structure factor at the 
lattice pinning point of the ground state = 0) for each 
lattice strength. 



and reaches constant when all vortices are pinned for 
Vo/hwo > 0.4. Correspondingly the structure factor 
|5'(K)| increases gradually and finally becomes constant 
when all vortices are pinned. From these results, together 
with directly looking at the condensate density profile, 
we conclude that the structural phase transition occurs 
at Vo/hooo = 0.4, which we define as the strength for the 
structural phase transition Vc/huJo. Fig. 3 (b) shows the 
analogous result for ft/ujo = 0.55. In this case, there is 
an intermediate domain where the Abrikosov lattice and 
the pinned lattice coexist. In this domain, the vortex lat- 
tice structure cannot be categorically determined because 
of the competition between the vortex- vortex interaction 
and the lattice pinning effect. 

From these results, we map out the phase diagrams of 
vortex lattice structures against ft and Vb, as we plot in 
Fig. 4 (a). The lower domain is the Abrikosov lattice 
domain, while the upper domain is the pinned lattice 
domain. The intermediate domain represents coexist- 
ing state of the Abrikosov lattice and the pinned lattice. 
Looking at the phase boundary of the pinned lattice do- 
main in Fig. 4 (a) as = K(r^), we find that Vc takes a 
minimum value as a function of which we define as the 
minimum pinning strength K,min = Vc{^ = ^^min)- From 
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FIG. 4: (a) Phase diagrams for vortex lattice structures for 
a/ttho = 2.2. (b) Lattice constant a against Qmin giving the 
minimum pinning strength. The solid line represents Eq. (5). 



Fig. 4 (a) for a/aho = 2.2, we find that K,min/^^o ~ 0.4 
and l^min/^o ~ 0.70. 

The dependence of minimum pinning lattice strength 
K,min on the lattice constant a can be understood as 
follows. When vortices form the triangular lattice and 
undergo rigid rotation with angular velocity 1], the lattice 
constant is expressed as a function of angular velocity by 



27r 



V3 2Vt/uoQ ' 



(5) 



In Fig. 4 (b), we plot lattice constant a against l^min- 
We find that the formula in Eq. ([5]) well fits the numeri- 
cal data. This means that the pinning strength Vc takes 
minimum value when a^(yt) matches to the lattice con- 
stant a. When this "matching relation" is satisfied, weak 
lattice potential has a stronger pinning effect than vortex- 
vortex interaction, which leads to a sharp transition of 
the vortex lattice structure. 

In the case when this "matching relation" is satisfied, 
one may presume that the vortex lattice array match 
the lattice potential array without vortex pinning effect. 
However, from Fig. 4 (a), we find that minimum pinning 
lattice strength K,min does not fall to zero, but takes a 
finite value. Actually in obtaining the phase diagram in 
Fiq. 4 (a), we solved the Gross-Pitaevski equation start- 
ing with the initial ground state wavefunction without 
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FIG. 5: (a) Phase diagrams for vortex lattice structures for 
a/aho = 2.2 by increasing Vb or decreasing Vo from pinned 
lattice domain (♦). At 11 > r^min, phase boundary of actually 
the pinned lattice phase is solid line, (b) Total free energy F 
with a/aho = 2.2 at Q/ooo = 0.70 (0), ^/(^o = 0.71 (Q) and 
Q/ooo = 0.725 (a). Here ^increase is in the case of increase of 
Vb , while ^decrease IS lu the casc of decrease of Vb from pinned 
lattice domain, (c) Total free energy F with a/aho = 2.2 at 
Q/loo = 0.65 (0), ^/(^o = 0.675 (Q) and Q/uo = 0.69 (a). 



vortices. As evolving the Gross-Pitaevski equation, vor- 
tices are nucleated, forming vortex lattices. In this dy- 
namical process, it may be possible that vortices relax 
into met ast able configurations, and this may be the rea- 
son why K,min takes a finite value. In order to investigate 
this in more detail, we solve the Gross-Pitaevski equation 
starting with a perfectly pinned vortex lattice at large 
enough Vq. We then slowly decrease Vq. In Fig. 5 (a). 



we plot the phase boundary of the pinned lattice phase 
obtained in this way, together with phase boundaries pre- 
viously shown in Fiq. 4 (a). We find that at 17 = l^min, 
the perfectly pinned vortex lattices is stable down to in- 
finitesimally small lattice potential Vb ^ 0. In order to 
see whether the pinned state obtained here has the lower 
energy than the tilted Abrikosov lattice obtained in the 
previous calculation, we calculate the total free energy 



F = 



|^V^ + Kxt(r)+^|V^(r,t)r 



(6) 



In Figs. 5 (b) and (c), we plot the free energy F for 
a/dho = 2.2 against Vq. Here, ^increase is the free energy 
obtained by solving the Gross-Pitaevski equation start- 
ing with non- vortex initial state and increasing Vq^ while 
^decrease IS the free energy obtained by decreasing Vb from 
pinned lattice domain. We make the separate plots for 
Vt > r^min and O < r^min in Figs. 5 (b) and (c), respec- 
tively. One can see from Fig. 5 (b) that, at = ^Imin, 

i^decrease < i^increase holds doWU tO ^ 0. This mcaUS 

that the Abrikosov lattice found for < < K,min in 
the phase diagram in Fig. 4 (a) is actually metastable. 
Similarly for Q > l^min, we find a finite region where the 
Abrikosov lattice is metastable. However, for Q < flmin^ 
we find ^decrease > ^increase as showu in Fiq. 5 (c). This 
means that the pinned vortex lattice does not have the 
lowest energy in the weak lattice domain. 

From the above results, we found that property of vor- 
tex pinning is quite different for Q > ftmin and Q < l^min- 
For Q > r^min, the vortex lattice structure exhibits 
sharp transition. In contrast, for Q < l^min, there is no 
sharp transition, but the vortex lattice structure exhibits 
crossover through the intermediate coexisting region. 



IV. SUMMARY AND DISCUSSION 

In summary, we have studied the vortex pinning effect 
by observing the structural phase transition of vortex lat- 
tice structures of a Bose-Einstein condensate in a rotating 
triangular lattice potential. We observed the transition 
of the vortex lattice structure from the Abrikosov vortex 
lattice to the pinned lattice where all vortices are pinned 
to lattice points. The transition is determined depending 
on the competition between strength and density of the 
lattice potential, and vortex density and vortex-vortex 
interaction. From our simulation, we found that the lat- 
tice potential has a strong vortex pinning effect when 
ay{Q) matches to the lattice constant a of the triangu- 
lar lattice potential, which means that the vortex den- 
sity coincides the density of pinning points. In the case 
when this "matching relation" is satisfied, we observed 
that minimum pinning lattice strength K,min takes a fi- 
nite value. By investigating the structural transition in 
more detail, we found that for Q > (^min, there are re- 
gions where the Abrikosov lattice is metastable, while for 
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Q < llmin, the pinned vortex lattice phase is metastable 
state. From the above results, we found that property 
of vortex pinning is quite different for ft > Omin and 
Q < Clmin- For Q > r^min, the vortex lattice structure ex- 
hibits sharp transition. In contrast, for Q < l^min, there 
is no sharp transition, but the vortex lattice structure ex- 
hibits crossover through the intermediate coexisting re- 
gion. 

In a separate paper, we will discuss the pinning effect 
of a rotating square lattice potential. We will show that 
the property of vortex pinning is quite different from the 
triangular lattice and more complicated for a rotating 



square lattice potential. 
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